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Abstract 

We analyse the phenomenology of orbifold scenarios from the heterotic su- 
perstring, and the resulting theoretical predictions for the direct detection of 
neutralino dark matter. In particular, we study the parameter space of these 
constructions, computing the low-energy spectrum and taking into account the 
most recent experimental and astrophysical constraints, as well as imposing the 
absence of dangerous charge and colour breaking minima. In the remaining al- 
lowed regions the spin-independent part of the neutralino-proton cross section is 
calculated and compared with the sensitivity of dark matter detectors. In addi- 
tion to the usual non universalities of the soft terms in orbifold scenarios due to 
the modular weight dependence, we also consider D-term contributions to scalar 
masses. These are generated by the presence of an anomalous f7(l), providing 
more flexibility in the resulting soft terms, and are crucial in order to avoid 
charge and colour breaking minima. Thanks to the D-term contribution, large 
neutralino detection cross sections can be found, within the reach of projected 
dark matter detectors. 
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1 Introduction 



One of the most interesting candidates for the dark matter in the Universe is a Weakly 
Interacting Massive Particle (WIMP), and in fact many underground experiments are 
being carried out around the world in order to detect its flux on the Earth [1]. These 
try to observe the elastic scattering of WIMPs on target nuclei through nuclear recoils. 
Although one of the experiments, the DAMA collaboration [2], reported data favouring 
the existence of a signal with WIMP-proton cross section 10~^ — 10~^ pb for a WIMP 
mass smaller than 500 — 900 GeV [2,3], other collaborations such as CDMS Soudan 
[4], EDELWEISS [5], and ZEPLIN I [6] claim to have excluded important regions of 
the DAMA parameter space^. Recently, the XENONIO experiment at the Gran Sasso 
National Laboratory [8] has set the strongest upper hmit for the WIMP-proton cross 
section, further disfavouring the DAMA result. This controversy will be solved in the 
future since many experiments are running or in preparation around the world. For 
example, LIBRA [9] and ANAIS [10] wiU probe the region compatible with DAMA 
result. Moreover, CDMS Soudan will be able to explore a WIMP-proton cross section 
(7 > 2 X 10~^ pb, and planned 1 tonne Ge/Xe detectors are expected to reach cross 
sections as low as 10~^° pb [11]. 

The leading candidate within the class of WIMPs is the lightest neutrahno, a 
particle predicted by supersymmetric (SUSY) extensions of the standard model. Given 
the experimental situation, and assuming that the dark matter is a neutralino, it is 
natural to wonder how big the cross section for its direct detection can be. This 
analysis is crucial in order to know the possibility of detecting dark matter in the 
experiments. In fact, the analysis of the neutralino-proton cross section has been car- 
ried out by many authors and during many years [1]. The most recent studies take 
into account the present experimental and astrophysical constraints on the parameter 
space. Concerning the former, the lower bound on the Higgs mass, the b —>■ sj and 
fi~^fi~ branching ratios, and the muon anomalous magnetic moment have been 
considered. The astrophysical bounds on the dark matter density, 0.1 < ^h"^ < 0.3 
(0.095 < Qh'^ < 0.112 if we take into account the recent data obtained by the WMAP 
satellite [12]), have also been imposed on the theoretical computation of the relic neu- 
tralino density, assuming thermal production. In addition, the constraints that the 
absence of dangerous charge and colour breaking minima imposes on the parameter 
space have also been implemented [13]. 

^For attempts to show that DAMA and these experiments might not be in conflict, see Ref. [7]. 
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In the usual minimal supergravity (mSUGRA) scenario, where the soft terms of 
the minimal supersymmetric standard model (MSSM) are assumed to be universal 
at the unification scale, Mqut ~ 2 x lO"*^^ GeV, and radiative electroweak symmetry 
breaking is imposed, the neutralino-proton cross section turns out to be constrained by 
cr^o_p < 3 X 10"^ pb. Clearly, in this case, present experiments are not sufficient and 
only the planned 1 tonne Ge/Xe detectors would be able to test part of the parameter 
space. However, in the presence of non-universal soft scalar and gaugino masses [14] 
the cross section can be increased significantly [15] in some regions with respect to the 
universal scenario (see, e.g., the discussion in [16], and references therein). Although 
the current upper limit on the decay Bg — > seriously affects these results, as 

was pointed out in [17], regions of the parameter space can still be found where the 
neutralino detection cross section can be within the reach of experiments such as CDMS 
Soudan. An analysis, summarizing all these results in the context of SUGRA, can be 
found in [18]. 

On the other hand, the low-energy limit of superstring theory is SUGRA, and 
therefore the neutralino is also a candidate for dark matter in superstring constructions. 
Let us recall that, in the late eighties, working in the context of the x heterotic 
superstring, a number of interesting four-dimensional vacua with particle content not 
far from that of the SUSY standard model were found (see, e.g., the discussion in 
the introduction of [19], and references therein). Such constructions have a natural 
hidden sector built-in: the complex dilaton field S arising from the gravitational sector 
of the theory, and the complex moduli fields Tj parametrizing the size and shape of 
the compactified space. The auxihary fields of those gauge singlets can be the seed 
of SUSY breaking, solving the arbitrariness of SUGRA where the hidden sector is not 
constrained. In addition, in superstrings the gauge kinetic function, fa{S,Ti), and the 
Kahler potential, K{S, S* ,Ti,T*), can be computed exphcitly, leading to interesting 
predictions for the soft parameters [20] . More specifically, in orbifold constructions they 
show a lack of universality due to the modular weight dependence. From these resulting 
SUGRA models one can also obtain predictions for the value of the neutralino-proton 
cross section. In fact, analyses of the detection cross section in these constructions 
were carried out in the past in [21-23]. 

Our aim in this work is to study in detail the phenomenology of these orbifold 
models, including the most recent experimental constraints on low-energy observables, 
as well as those coming from charge and colour breaking minima, and to determine 
how large the cross section for the direct detection of neutralino dark matter can be. 
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We therefore calculate the theoretical predictions for the spin-independent part of the 
neutralino-nucleon cross section, O'^o_p, and compare it with the sensitivities of present 
and projected experiments. Since the soft terms in superstring scenarios are a subset 
of the general soft terms studied in SUGRA theories we make use of previous results 
on departures from the mSUGRA scenario to look for values of the orbifold soft terms 
giving rise to a large cross section accessible for experiments. 

In addition, we introduce a new ingredient in the analysis, namely the modification 
produced in the soft parameters by the presence of an anomalous U{1). Let us recall 
that in string theory, and in particular in orbifold constructions [24, 25] of the heterotic 
superstring [26], the gauge groups obtained after compactification are larger than the 
standard model gauge group, and contain generically extra U{1) symmetries, SU{3) x 
SU (2) X U (1)" [27]. One of these t/(l)'s is usually anomalous, and although its anomaly 
is cancelled by the four-dimensional Green-Schwarz (GS) mechanism, it generates a 
Fayet-IHopoulos (FI) contribution to the D-term [28]. This effect is crucial for model 
building [29] since some scalars acquire large vacuum expectation values (VEVs) in 
order to cancel the FI contribution, thereby breaking the extra gauge symmetries, and 
allowing the construction of realistic standard-like models in the context of the 
orbifold [30-32] (see also [33]). Recently other interesting models in the context of 
the Ze orbifold [34-36], and Z12 orbifold [37,38], have been analysed. Due to the FI 
breaking, also D-term contributions to the soft scalar masses are generated [39-44]. 
This allows more flexibility in the soft terms and, consequently, in the computation of 
the associated neutralino-proton cross section. 

The paper is organised as follows. In Section 2 we briefly review the departures 
from mSUGRA which give rise to large values of the neutralino detection cross section. 
Then, in the next sections, we use this analysis to study several orbifold scenarios 
where such departures may be present. Special emphasis is put on the effect of the 
various experimental constraints on the SUSY spectrum and low-energy observables. 
We start in Section 3 with the simplest (but not the most common) possibility, where 
an anomalous U{1) is not present. Then, in Section 4, we discuss the important 
modifications produced in the soft terms by the presence of an anomalous U{1), and 
their effects on the computation of the neutralino-proton cross section, considering 
the effect of D-term contributions to soft scalar masses. The conclusions are left for 
Section 5. 
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2 Neutralino-proton cross section and departures 
from mSUGRA 

In this section we review possible departures from the mSUGRA scenario, and their 
impact on the neutrahno-proton cross section. This will allow us to discuss orbifold 
scenarios more easily. Let us first recall that in mSUGRA one has only four free pa- 
rameters defined at the GUT scale: the soft scalar mass, m, the soft gaugino mass, M, 
the soft trilinear coupling, A, and the ratio of the Higgs vacuum expectation values, 
tan/5 = {H^)/{H^). In addition, the sign of the Higgsino mass parameter, fi, remains 
undetermined. Using these inputs the neutralino-proton cross section has been anal- 
ysed exhaustively in the literature, as mentioned in the Introduction. Taking into 
account all kind of experimental and astrophysical constraints, the resulting scalar 
cross section is bounded to be cr^o_p < 3 x 10~® pb. 

Departures from the universal structure of the soft parameters in mSUGRA allow 
to increase the neutralino-proton cross section significantly. As it was shown in the 
literature, it is possible to enhance the scattering channels involving exchange of CP- 
even neutral Higgses by reducing the Higgs masses, and also by increasing the Higgsino 
components of the lightest neutralino. A brief analysis based on the Higgs mass pa- 
rameters, rn'jj^ and m%^, at the electroweak scale can clearly show how these effects 
can be achieved. 

First, a decrease in the values of the Higgs masses can be obtained by increasing 
in%^ at the electroweak scale (i.e., making it less negative) and/or decreasing mjj^. 
More specifically, the value of the mass of the heaviest CP-even Higgs, H, can be very 
efficiently lowered under these circumstances. This is easily understood by analysing 
the (tree-level) mass of the CP-odd Higgs A, which for reasonably large values of tan f3 
can be approximated as m\ ^ mf^^ —m^^ — M|. Since the heaviest CP-even Higgs, H , 
is almost degenerate in mass with A, lowering m\ we obtain a decrease in m\ which 
leads to an increase in the scattering channels through Higgs exchange 

Second, through the increase in the value of rn^j^ an increase in the Higgsino com- 
ponents of the lightest neutralino can also be achieved. Making less negative, its 
positive contribution to /x^ in the minimization of the Higgs potential would be smaller. 
Eventually will be of the order of Mi, M2 and x5 ^i^^ then be a mixed Higgsino- 
gaugino state. Thus scattering channels through Higgs exchange become more impor- 
tant than in mSUGRA, where \^\ is large and x5 is mainly bino. It is worth emphasizing 
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however that the effect of lowering the Higgs masses is typically more important, since 
it can provide large values for the neutralino-nucleon cross section even in the case of 
bino-like neutralinos. 

Non-universal soft parameters can produce the above mentioned effects. Let us 
consider in particular the non-universality in the scalar masses, which will be the most 
interesting possibility in orbifold scenarios. We can parametrize these in the Higgs 
sector, at the high-energy scale, as follows: 

^ m\l + S„,) , ml^^m'{l + 5Hj . (2.1) 

Concerning squarks and sleptons we will assume that the three generations have the 
same mass structure: 

= m\l + 6Lj. (2.2) 

Such a structure avoids potential problems with flavour changing neutral currents^ 
(FCNC), and arises naturally e.g. in orbifold constructions with two Wilson lines, 
where reahstic models have been obtained. Note also that whereas all 5's in (2.2) have 
to satisfy 5 > — 1 in order to avoid an unbounded from below (UFB) direction breaking 
charge and colour^, Sh^^ < —1 in (2.1) is possible as long as ml — mjj^ + /j,'^ > and 
m| = rriH^ + /j,'^ > are fulfilled. 

An increase in m|^^ at the electroweak scale can be obviously achieved by increasing 
its value at the high-energy scale, i.e., with the choice Sh„ > 0. In addition, this is also 
produced when rriq^ and m^^ at the high-energy scale decrease, i.e. taking 5q^^ur < 0, 
due to their (negative) contribution proportional to the top Yukawa coupling in the 
renormalization group equation (RGE) of rriff^. 

Similarly, a decrease in the value of m|^^ at the electroweak scale can be obtained by 
decreasing it at the high-energy scale with Sh^ < 0. The same effect is obtained when 
rriQ^ and m^^ increase at the high-energy scale, due to their (negative) contribution 

^Another possibility would be to assume that the first two generations have the common scalar 
mass m, and that non-universalities are allowed only for the third generation (as it occurs for the 
models analysed in Ref. [45]). This would not modify our analysis since, as we will see below, only the 
third generation is relevant in our discussion. 

■^If we allow metastability of our vacuum, tachyonic masses for some sfcrmions, S < —1, at the 
high-energy scale might be allowed. However, we do not consider such a possibility. 
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proportional to the bottom Yukawa coupling in the RGE of Thus one can deduce 
that m\ will also be reduced by choosing Sq^^Hr > 0. 

In fact non-universality in the Higgs sector gives the most important effect, and 
including the one in the sfermion sector the cross section only increases slightly. 

Taking into account this analysis, several scenarios were discussed in Ref. [13], 
obtaining that large values for the cross section are possible. For example, with Sh^ — 
0, = 1; — —1, = 0; — —1, ^Hu — 1) regions of the parameter space are 
found which are accessible for experiments such as CDMS Soudan [18]. Interestingly, 
it was also realised that these choices of parameters were helpful in order to prevent 
the appearance of UFB minima in the Higgs potential. 

The different UFB directions were classified in Ref. [46]. Among these, the one 
labelled as UFB-3, which involves VEVs for the fields {Hu^VLiT^Lj^^Rj} with i ^ j, 
yields the strongest bound. After an analytical minimization of the relevant terms of 
the scalar potential the value of the VLi, cl^, ^.r. fields can be written in terms of Hu- 
Then, for any value of \Hy\ < Mqut satisfying 



4X1 



+ 



-,I2 



9l 



2A, 



(2.3) 



the potential along the UFB-3 direction reads 



Otherwise 



Kjfb-3 = <Ji^«P + ^-^{m% + ml)\H^\ + \{g'' + gl) 



2mi. 



1 2 



\Hu\ + —\Hu\ 



(2.4) 



(2.5) 



In these expressions Ae^ denotes the leptonic Yukawa coupling of the jth generation, the 
deepest direction corresponding to ej — r. The UFB-3 condition is then VupBsiQ = 

Q) > Keal min.(Q = MsuSy), where Vreal min. = -| (g'"^ + gl) {vl - V^f , with Vu,d = 

{Hu,d), is the value of the potential at the realistic minimum. Vreai min is evaluated at 
the typical scale of SUSY masses, Msusy, and Vufb-s at the renormahzation scale, 
(5, which is chosen to be Q ~ Max(Atop|Hu|, Mgusv), in order to minimize the effect of 
one-loop corrections to the scalar potential. 

As we see from Eqs. (2.4) and (2.5), the potential along this direction can be lifted 
when m\j increases (becomes less negative) and for large values of the stau mass 
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parameters, thereby making the UFB-3 condition less restrictive. In this sense, non- 
universal soft terms, like the ones discussed above, can be very helpful. 

The question now is whether it is possible to find explicit realisations of these 
scenarios within orbifold models. In the following sections we will study this issue in 
detail. 

3 Orbifold scenarios 

Let us recall first the structure of the SUGRA theory in four-dimensional constructions 
from the heterotic superstring. The tree-level gauge kinetic function is independent of 
the moduli sector and is simply given by 

fa = kaS , (3.6) 

where ka is the Kac-Moody level of the gauge factor. Usually (level one case) one 
takes ks — k2 — f/ci = 1 for the MSSM. In any case, the values ka are irrelevant for 
the tree-level computation since they do not contribute to the soft parameters. On 
the other hand, the Kahler potential has been computed for six-dimensional Abelian 
orbifolds, where three moduli Tj are generically present. For this class of models the 
Kahler potential has the form 

K = -log{S + S*)-J2^og{Ti + T*) + J2\C''\''ni{Ti + T:)< . (3.7) 

i a 

Here are (zero or negative) fractional numbers usually called 'modular weights' of 
the matter fields C". 

In order to determine the pattern of soft parameters it is crucial to know which 
fields, either S or Ti, play the predominant role in the process of SUSY breaking. Thus 
one can introduce a parametrization for the VEVs of dilaton and moduli auxiliary fields 
[47]. A convenient one is given by [47-49] 

F' = V3{Ti + T*)m3/2COseei , (3.8) 

where i = 1,2,3 labels the three complex compact dimensions, 1713/2 is the gravitino 
mass, and the angles 9 and 0j, with |0.tp = 1, parametrize the Goldstino direction 
in the S, Ti field space. Here we are neglecting phases and the cosmological constant 
vanishes by construction. 
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Using this parametrization and Eqs. (3.6) and (3.7) one obtains the following results 
for the soft terms [47-49]: 

Ma = 777.3/2 sin 6* , 

+ coseJ2^i[^+ < + ni3 + ni,-{Ti + T*)dilogX^p^]j . (3.9) 

Although in the case of the A parameter an explicit Tj-dependence may appear in 
the term proportional to dilogXaf^-y, where Xa/S'yiTi) are the Yukawa couplings and 
di = d/dTi, it disappears in several interesting cases [47,49]. For example, the A- 
term which is relevant to electroweak symmetry-breaking is the one associated to the 
top-quark Yukawa coupling. Thus, in order to obtain the largest possible value of the 
coupling, the fields should be untwisted or twisted associated to the same fixed point. 
In both cases diXa^e-y 0, and we will only consider this possibility here. 

Using the above information, one can analyse the structure of soft parameters 
available in Abelian orbifolds. In the dilaton-dominated SUSY-breaking case (cos^ = 
0) the soft parameters are universal, and fulfil [50,51] 

777 = 7713/2 , M = ±v^m , A = -M , (3.10) 

where the positive (negative) sign for M corresponds to 9 — n /2 (9 — Sn /2) . Of course, 
these are a subset of the parameter space of mSUGRA, and as a consequence one should 
expect small dark matter detection cross sections, as discussed in the previous Section. 

However, in general, the soft terms (scalar masses and trilinear parameters) given 
in Eq. (3.9) show a lack of universality due to the modular weight dependence. For 
example, assuming an overall modulus (i.e., T — and ©j = l/y/3), one obtains 

"^a = 7773/2 (l +7l„COS^ 6') , (3.11) 
Aap-y = -"\/3 7773/2 sin ^ - 777.3/2 COS ^ (3 7la 71^ 71-y) , (3-12) 

where wc have defined the overall modular weights Ua = Yli'^^a- ^^e case of Z„ 
Abelian orbifolds, these can take the values —1, —2, —3, —4, —5. Fields belonging to 
the untwisted sector of the orbifold have Ua — —1. Fields in the twisted sector but 
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without oscillators have usually modular weight —2, and those with oscillators have 
iT'a ^ ~3. Of course, if all modular weights of the standard model fields are equal, 
one recovers the universal scenario. For example, taking all na = — 1 one has [47] 
m = ms/2sm9, M = \^m, A = —M. 

Using notation (2.1) and (2.2), the degree of non-universality in the scalar masses 
is therefore given by 

Sa = ria cos^ . (3.13) 

It is worth noticing here that 5q, < as a consequence of the negativeness of the 
modular weights. As we will see, this has important phenomenological implications. 

On the other hand, the apparent success of the joining of gauge coupling constants 
at, approximately, 2 x 10^® GeV in the MSSM is not automatic in the heterotic super- 
string, where the natural unification scale is Mqut — gauT x 5.27 x 10^^ GeV, where 
qgut is the unified gauge coupling. Therefore unification takes place at energies around 
a factor 10 smaller than expected in the heterotic superstring. This problem might be 
solved with the presence of large string threshold corrections which explain the mis- 
match between both scales [52,53]. In a sense, what would happen is that the gauge 
coupling constants cross at the MSSM unification scale and diverge towards different 
values at the heterotic string unification scale. These different values appear due to 
large one-loop string threshold corrections. 

It was found that these corrections can be obtained for restricted values of the 
modular weights of the fields [53]. In fact, assuming generation independence for the 
ria as well as —3 < n„ < —1, the simplest possibility corresponds to taking the following 
values for the standard model fields: 

'^Ql = "^dR = -1, nuji = -2, uli^ = rie^ = -3 , 

^/f. + rij/, = -5, -4 , (3.14) 

where, e.g., u^i denotes the three family squarks -U/?, c/j, t^. The above values together 
with ReT ~ 16 lead to good agreement for sin^ 6'vk and [53]. The associated soft 
sfermion masses are given by [47]: 

2 2 2 /1 2 a\ 

^Ql^ ^dR = 7723/2(1 - COS 9) , 

^Ir = ^3/2 (1-2 cos^ 9) , 
^I^^^Ir = m^/2(l-3cos2^) , (3.15) 
whereas for the soft Higgs masses, choosing uh^ = — 1 , uh^ — —3, one obtains: 



10 





















1 

- i 


9 


1 

- i 


-o 


Q 

-o 


Q 

-o 


-i 


B) 


-1 


-2 


-1 


-3 


-3 


-3 


-2 


C) 


-1 


-2 


-2 


-1 


-1 


-2 


-1 


D) 


-2 


-1 


-1 


-2 


-1 


-2 


-2 


E) 


-1 


-1, -3, -3 


-1 


-3 


-1, -3, -3 


-2 


-3 



Table 1: Modular weights for the scalar fields of heterotic orbifold scenarios with an overall 
modulus that can reproduce gauge unification [50,53]. Note that cases B) and C) present 
extra massless chiral fields. 

niH^ = ml/2 (1-3 cos^ 9) . (3.16) 

For convenience, this set of modular weights is summarised in Table 1 and labelled as 
case A). 

For example, with cos^ ^ = 1/3, using notation (2.1) and (2.2), the non-universalities 
in the Higgs and sfermion sectors correspond to Sh^, — —1/3, Sh^ — —1; ^Ql — — 
-1/3, 5uji = -2/3, and 5lj^ = = -1- 

Concerning the soft gaugino masses, they are given by: 

Ms ~ l.oVSmg/ssin^ , 

Ml ~ 1.18\/3m3/2sin^ . (3.17) 

The small departure from universality is due to the effect of the string threshold cor- 
rections on the gauge kinetic function [47] . 

Finally, for the above modular weights, and using (3.12), the expressions for the 
trilinear parameters read 

Aj- — — m3/2(\/3 sin6' — 6 cos6') , 
A}, — — m3/2(-\/3sin^ — 2 cos^) , 

At = -m3/2(V3sin^ -cos^) . (3.18) 

The yl-term which is relevant to radiative symmetry breaking is the one associated to 
the top-quark Yukawa coupling At. 

These soft terms serve as an explicit model for the study of the ncutralino detection 
cross section. Since they are completely determined in terms of just the gravitino mass 
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7t/2 71 37t/2 27T 

Figure 1: Soft terms at the string scale in units of 7713/2 as a function of the Goldstino 
angle, 6. Solid lines represent, from top to bottom, the bino, wino and gluino mass 
parameters. The various scalar masses are depicted by means of dashed lines. Finally, 
dotted lines correspond to the trilinear terms. The oblique ruled areas are excluded due to 
the occurrence of negative mass-squared parameters. 

and the Goldstino angle, we are left with three free parameters, namely m3/2, d, and 
tan/3, plus the sign of fi. Note, however, that the absence of negative mass-squared 
of the sleptons at the GUT scale implies the constraint cos^ < |. Besides, the shift 
6 ^ 6 + 71 implies in the above equations nia rria, Ma —Ma and At —At. This 
fact makes it unnecessary to consider both signs of the /i parameter. The reason is that 
the RGEs arc symmetric under the change /i, M, A — > — /i, — M, —A. Consequently, 
in the remainder of this paper we will assume /i > 0. Notice in this sense that we will 
always have yuMj > for < vr whereas /xMj < for 6 > n. This will have important 
implications, as we will soon see, on the effect of the experimental constraints on the 
rare decays 6 — > 57 and 5° f^^l^~, and on the SUSY contribution to the muon 
anomalous magnetic moment, a^^^ . 

The resulting structure of the soft parameters for case A), given at the GUT scale, 
is represented in Fig. 1 as a function of the Goldstino angle in units of the gravitino 
mass. Two generic features of this kind of orbifold constructions are evidenced by the 
plot, namely, the fact that scalar masses arc always smaller than gaugino masses, and 
the presence of regions which are excluded because some scalar masses-squared become 
negative. In the present example, as already mentioned, the strongest bound is set by 
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Figure 2: Supersymmetric spectrum at low-energy as a function of the Goldstino angle, 
6, for 777.3/2 = 200 GeV and tan/3 = 10 and 35. Only the region around 9 — 7r/2 is 
represented. From bottom to top, the solid lines correspond to the lightest neutralino, 
lightest chargino and gluino masses. The dashed lines represent the lightest stau and 
lightest sneutrino masses. The lightest stop and sbottom masses are plotted by means 
of dot-dashed lines. Finally, the dotted lines show the mass of the CP-odd Higgs and 
the resulting fj, parameter. The oblique ruled areas are excluded due to the occurrence of 
negative mass-squared parameters at the GUT scale, whereas the gridded regions correspond 
to those where tachyons appear after solving the RGBs. 

slepton masses, for which (3.15) imphes cos^ ^ < 1/3. The ruled areas correspond to 
those where this bound is not fulfilled. This reduces the allowed parameter space to 
two strips in 9, around the dilaton-dominated case, ^ = 7r/2, 37r/2. 

With this information, the RGEs are numerically solved and the low-energy super- 
symmetric spectrum is calculated. Fig. 2 shows the resulting particle spectrum as a 
function of the Goldstino angle for 771,3/2 = 200 GeV and tan/3 = 10 and 35. As we can 
see, although slepton masses-squared are positive at the GUT scale for cos^^ < 1/3, 
the RGEs can still drive them negative, or lead to tachyonic mass eigenstates. This is 
typically the case of the lightest stau, fi, and lightest sneutrino, i>i (the latter only for 
low values of the gravitino mass), due to their small mass parameters (3.15). This is 
more likely to happen for large tan f3, since the lepton Yukawas (which are proportional 
to 1/ cos /3) increase and induce a larger negative contribution to the slepton RGEs. In 
such a case, the lightest stau can be the lightest SUSY particle (LSP) in larger regions 
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of the parameter space, thus potentially reducing the allowed areas for neutralino dark 
matter, as we see in the example with tan/? = 35. The supersymmetric spectrum 
also displays a heavy squark sector, due to the gluino contribution on the running of 
their mass parameters. Similarly, the heavy Higgs masses (represented here only with 
the pseudoscalar, ^4°) are also sizable. For reference, the value of the /j, term is also 
displayed and found to be large. 

Notice at this point that there are regions of the parameter space where the lightest 
neutralino is the LSP and the stau, being the next-to-lightest SUSY particle (NLSP), 
has a very similar mass. As we will soon see, this allows reproducing the correct dark 
matter relic density by means of a coannihilation effect. On the other hand, one can 
readily see that in these examples m^o > 2 m^o and therefore there is no enhancement 
in the annihilation of ncutralinos mediated by the CP-odd Higgs. 

Finally, it is worth emphasizing that in these scenarios the gravitino is never the 
LSP. Despite the bino mass being larger than 7713/2 at the string scale, its RGE always 
leads to Mi < 771^/2 at the electroweak scale (even in the dilaton-dominated limit for 
which Ml is at its maximum) so that, at least, the neutralino mass is always lighter 
than 777.3/2. 

Having extracted the supersymmetric spectrum, we are ready to determine the 
implications for low-energy observables and study how the associated bounds further 
restrict the allowed parameter space. In our analysis the most recent experimental 
and astrophysical constraints have been included. In particular, the lower bounds 
on the masses of the supersymmetric particles and on the lightest Higgs have been 
implemented, as well as the experimental bound on the branching ratio of the b — > 
S7 process, 2.85 x 10^^ <B(6 s'~f)< 4.25 x 10^"^. The latter has been calculated 
taking into account the most recent experimental world average for the branching 
ratio reported by the Heavy Flavour Averaging Group [54-56], as well as the new re- 
evaluation of the SM value [57], with errors combined in quadrature. We also take into 
account the improved experimental constraint on the n^fj,^ branching ratio, 

B(i?° l^^l^ ) < 1-5 X 10~^, obtained from a combination of the results of CDF 
and DO, [58-60]. The evaluation of the neutralino rehc density is carried out with the 
program micrOMEGAs [61], and, due to its relevance, the effect of the WMAP constraint 
will be shown explicitly. Finally, dangerous charge and colour breaking minima of the 
Higgs potential will be avoided by excluding unbounded from below directions. 

Concerning a^^^^, we have taken into account the experimental result for the muon 
anomalous magnetic moment [62], as well as the most recent theoretical evaluations of 
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the Standard Model contributions [63-65]. It is found that when e^e~ data are used 
the experimental excess in = ((^^ — 2)/2 would constrain a possible supersymmetric 
contribution to be a^^^^ = (27.6 ± 8) x 10^^°, where theoretical and experimental 
errors have been combined in quadrature. However, when tau data arc used, a smaller 
discrepancy with the experimental measurement is found. Due to this reason, in our 
analysis we will not impose this constraint, but only indicate the regions compatible 
with it at the 2a level, this is, 11.6 x 10"^° < a^^SY < 43.5 x 10"^°. 

For a better understanding of all these constraints, we have represented in Fig. 3 
their effect on the (^3/2,^^) plane for tan/5 = 10. For comparison, the cases with 
tan f3 = 20, 35 are also shown in Fig. 4. 

The first thing to notice is that extensive regions are excluded due to the occurrence 
of tachyonic masses for sleptons. As already discussed, the area excluded for this reason 
becomes larger when tan/5 increases, an effect which is clearly displayed in Figs. 3 and 
4. This implies an increase in the lower bound of the gravitino mass. Whereas for 
tan (3 = 10, 20 the smallest allowed value is 777.3/2 ~ 35 GeV, in the case with tan /3 = 35 
one needs 1113/2 > 90 GeV. 

The above mentioned smallness of the slepton mass parameters, together with the 
fact that gaugino masses are always larger than scalar masses (M^ > rria), also imply 
that the areas in the parameter space where the lightest neutralino is the LSP are 
not very extensive. These regions occur for small values of cos^^ (they are centered 
around 9 ^ 7r/2, 37r/2), since the ratio \M\/mLj^^eR increases with"^ cos^^. Note that 
such values of the Goldstino angle mean that the breaking of SUSY is mainly due to 
the dilaton auxiliary term. Once more, the allowed areas shrink for large vahics of 
tan/3 and eventually disappear for tan/5 ^ 45. In the rest of the parameter space the 
role of the LSP is mainly played by the lightest stau. Although, as already mentioned, 
for small values of tan (3 the sneutrino can also be the LSP in a very narrow band for 
small gravitino masses, this area is always excluded by experimental bounds. 

The relevance of the experimental constraints is also evidenced by Figs. 3 and 4. 
Reproducing the experimental result of the branching ratio of 6 ^ 57 is much easier 
in the region around 6 = 7r/2, since it has fiM > 0. On the contrary, it poses a 
stringent lower bound on the value of 7/23/2 for the region around 9 = 37r/2, for which 
HM < 0. As expected, the area excluded for this reason also increases for larger values 
of tan/3. Thus, whereas this constraint implies 7/23/2 > 150 GeV for tan/5 = 10 in the 

"^The lack of a complete mirror symmetry at ^ = tt/2 and 37r/2 is due to the trilinear terms (3.12) 
being a combination of s'm9 and cos 9. 
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7t/4 7y/2 37t/4 57t/4 37t/2 77t/4 

Figure 3: Effect of the different experimental constraints on the parameter space {ms/2, 0) 
for tan/3 = 10. The oblique ruled areas are excluded due to the occurrence of tachyons 
at the GUT scale, whereas the gridded regions correspond to those where tachyons appear 
after solving the RGBs. Only the two areas centred around 9 ~ 7r/2, 37r/2 are free from 
tachyons in the slepton sector. Dark grey areas represent those where the lightest neutralino 
is not the LSP. Among these regions, the narrow vertical areas contained within solid lines 
are those where the stau is the LSP, whereas in the thin horizontal region at very small 
gravitino masses, also bounded by solid lines, the LSP is the lightest sneutrino. Light 
grey areas stand for those not fulfilling one or more experimental bounds. In particular, 
the region below the thin dashed line is excluded by the lower bound on the Higgs mass. 
The area below the thick dashed line is excluded by 6 57. The regions excluded by 
the experimental constraints on the masses of the chargino and stau are always contained 
within those ruled out by other reasons and are therefore not shown. The region bounded 
by thin dot-dashed lines is favoured by a^^^ (notice that the whole allowed area around 
9 ~ 3n/2 always has a^^^^ < and is therefore disfavoured), although this constraint has 
not been explicitly applied. The remaining white area is favoured by all the experimental 
constraints. Within it the ruled region fulfils in addition 0.1 < fl^oh'^ < 0.3, and the black 
area on top of this indicates the WMAP range 0.094 < Vl^oh^ < 0.112. Finally, the UFB 
constraints are only fulfilled in the area above the thick solid line. 

area around 9 — 7r/2, 777,3/2 > 250 (300) GeV is necessary for tan/3 = 20 (35). 

Having fiM < 0, the whole region around 9 = Sir/ 2 also fails to fulfil the experi- 
mental constraint on a^^^^, and is therefore further disfavoured. 
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1? 

Figure 4: The same as Fig. 3 but for tan/3 = 20 and 35. Notice that the UFB constraints 
are violated in the whole represented parameter space and therefore no thick solid line is 
shown. In the plot for tan/3 = 35, the area below the thick dot-dashed line on the region 
is excluded due to the constraint on the — > branching ratio 

The bound on the hghtest Higgs mass also rules out some regions for small gravitino 
masses. This is only relevant for small values of tan f3 and in the region around 9 = 7r/2. 
Already for tan /? > 15 this bound becomes less important than the 6 — > 57 or a^^^^ 
constraints. The areas not fulfilling the experimental constraints on sparticle masses 
are always contained within those already excluded by other bounds and are therefore 
not shown explicitly. 

The allowed parameter space is further reduced when the constraint on the relic 
density is imposed. The WMAP result is only reproduced along the narrow regions 
close to the area where the stau becomes the LSP. This is due to the well known 
coannihilation effect that takes place when the neutralino mass is close to the stau mass. 
The equivalent of the "bulk region" in the mSUGRA parameter space is here excluded 
by the experimental constraints. Finally, no regions are found where 2m^o ^ m^, and 
consequently resonant annihilation of neutralinos does not play any role in this case. 

Having shown that there are regions with viable neutralino dark matter, let us now 
turn our attention to its possible direct detection. Following the discussion of Sec- 
tion 2, the Higgs modular weights giving rise to the soft masses (3.16), could induce an 
increase of the neutralino detection cross section with respect to the universal case. In 
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Figure 5: Scatter plot of the scalar neutralino-proton cross section, o'^o_p, as a function 
of the neutralino mass, m^o, for tan/3 = 10, 20 and 35. The light grey dots represent 
points fulfilling all the experinnental constraints. The dark grey dots represent points which 
satisfy in addition 0.1 < fl^oh'^ < 0.3 and the black dots on top of these indicate those 
in agreement with the WMAP range 0.095 < ^^oh"^ < 0.112. The sensitivities of present 
and projected experiments are also depicted with solid and dashed lines, respectively, in the 
case of an isothermal spherical halo model. The large (small) area bounded by dotted lines 
is allowed by the DAMA experiment when astrophysical uncertainties to this simple model 
are (are not) taken into account. 

order to investigate this possibility, the theoretical predictions for the spin-independent 
part of the neutralino-nucleon cross section have been calculated in the accepted re- 
gions of the parameter space. They are represented versus the neutralino mass in 
Fig. 5 for tan/3 = 10, 20 and 35, where the sensitivities of present and projected dark 
matter experiments are also shown. These results resemble those of mSUGRA, as no 
high values are obtained. As in mSUGRA, in this scenario the fj, parameter and the 
heavy Higgs masses are sizable (see Fig. 2), thus implying bino-like neutralinos and a 
suppressed contribution to cr^o_p from Higgs-exchanging processes. This is illustrated 
in Fig. 6, where the resulting values of the /j, parameter are plotted as a function of the 
pseudoscalar Higgs mass. After analysing the range tan /? = 10 to 50 we found that 
o'xi-p ~ ^ ^ 10~^ pb, the maximum values corresponding to tan (3 ^ 15. These results 
are therefore beyond the present sensitivities of dark matter detectors and would only 
be partly within the reach of the projected 1 tonne detectors. 
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Figure 6: Scatter plot of the resulting ji parameter as a function of the pseudoscalar Higgs 
mass, rriA, for tan/3 = 10, 20 and 35. The colour convention is as in Fig. 5. 

So far we have not commented on the bounds imposed by the UFB-3 constraint to 
avoid dangerous charge and colour breaking minima of the Higgs potential. This turns 
out to play a crucial role in disfavouring this scenario. Indeed, most of the parameter 
space is excluded on these grounds^. Only for small values of tan (5 and heavy gravitinos 
do allowed regions appear (see for instance Fig. 3, where 7713/2 ^ 650 GeV is necessary), 
but these always correspond to areas where the neutralino relic density is too large and 
exceeds the WMAP constraint. For tan/3 > 15 the UFB-3 constraint already excludes 
the complete region with 777.3/2 < 1000 GeV. Once more, the reason for this is the low 
value of the slepton masses, and more specifically, of the stau mass. Let us recall that 
the smaller this value, the more negative Vufb-z in (2.4) or (2.5) is, and thus the 
stronger the UFB-3 bound becomes. Moreover, the fact that in this scenario the value 
of rn^Hu ^®^y large (since 5hu is negative) also contributes in driving the potential 

deeper along this direction. 

Let us finally remark that other examples with different choices of modular weights 
for the Higgs parameters satisfying (3.14) have been investigated, such as case B) in 
Table 1, and lead to quahtatively similar results. 

The previous analysis suggests how to modify the model to 'optimise' its behaviour 
under the UFB-3 constraint [66], increasing also the regions in the parameter space 

^This is consistent with previous analyses of charge and colour breaking minima in different super- 
string and M- theory scenarios [66] . 
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Figure 7: The same as Fig. 2 but for the optimised example where = = —1. 

where the hghtest neutrahno is the LSP. The most favorable case would correspond to 
slepton masses as large as possible, i.e., 

riLL = neii = -1 ■ (3-19) 

For squark and Higgs mass parameters we will continue using the modular weights of 
case A)^. Note that now the bound on cos^ 9 is less constraining, since we only need 
to impose cos^ ^ < |, thus allowing a larger degree of non- universality. For example, 
with cos^ ^ = |, we get Sh^ — —1/3 and Sh^ = — 1 for the Higgs masses. 

The resulting super symmetric spectrum is shown in Fig. 7 as a function of the 
Goldstino angle for 771,3/2 = 200 GeV with tan/3 = 10 and 35. Notice that now the 
whole region with 9 e [7r/4, 37r/4] is free from tachyons at the GUT scale. For small 
tan/3 the increase in the slepton mass-squared parameters leaves extensive allowed 
regions where the neutralino can be the LSP. As expected, larger values of tan (3 lead 
to a reduction in the stau mass, which now easily becomes the LSP, and gives rise to 
tachyons in some regions. 

These features are evidenced in Fig. 8, where the corresponding (777.3/2, 9) parameter 

^Of course, with such a choice of modular weights we know that the string threshold corrections 
cannot account for the joining of gauge couplings at the MSSM unification scale. Thus we will be 
tacitly assuming that there is some other effect (e.g., the existence of further chiral fields in the 
spectrum below the hctcrotic string scale [19,33,69]) which appropriately produces the correct low- 
energy experimental values for gauge couplings. 
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Figure 8: The same as Fig. 3 but for the optimised example where n^,^ = rie^ = —1 with 
tan/3 = 20, 35. 



space is depicted for tan (3 = 20, 35. Notice that, due to the increase in the slepton 
mass terms, the stau only becomes the LSP on narrow bands on the right-hand side of 
the allowed areas for tan ^ — 20 (for smaller values of tan (3 the neutralino is always the 
LSP) . As expected, these areas with stau LSP become more sizable for tan /3 = 35 and 
eventually dominate the whole parameter space for tan/? > 45. Also, for tan/? < 25 
the sncutrino can also be the LSP on a narrow region for very light gravitinos, although 
this is always excluded by experimental constraints. 

The decrease of the stau mass towards the right-hand side of the allowed areas can 
be understood by analysing the expressions for the trilinear soft terms. The trilinear 
terms associated to the top, bottom and tau Yukawa coupling read in this example 

At — — 1712/2 (Vs sin 9 — cos 9^ , 
Ab^ = —ms/2 (V3 sin ^ — 2 cos . (3.20) 

It can be checked that for all of them the ratio |A/M| increases towards the right- 
hand side of both allowed areas. In particular, \Arfi/M\ ^ 0.05 for 9 = 7r/4, 57r/4 
and becomes \AT^b/M\ ^ —1.8 for 9 = 37r/4, 77r/4. The increase in |A^^f,/M| leads to 
a larger negative correction in the RGE for the slepton mass terms, implying lighter 
staus. Large values of tan P increase the corresponding Yukawas thus further decreasing 
the stau mass. 
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The variation in the stau mass affects the area excluded by the UFB-3 constraint, 
which becomes more stringent towards the right-hand side of the allowed regions. On 
the left, as expected, the effect of the UFB constraints is less severe than in the previous 
examples and regions with 7713/2 > 250 (400) GeV for tan P — 20 (35) are allowed. 
Interestingly, for tan /3 > 30 part of these areas can also reproduce the correct value 
for the neutralino relic density. 

The corresponding predictions for the neutralino-nucleon cross section are depicted 
in Fig. 9. Although regions with the correct rehc density can appear with cr^o_p Ki 10~^ 
pb for tan /3 = 20, the points fulfilling the UFB-3 constraints only correspond to those 
with a^o_p < 10~^ pb for tan /3 > 30. Once more, as evidenced in Fig. 10, where the 
resulting /i parameter is represented versus the CP odd Higgs mass, this is due to the 
large values of and the heavy Higgs masses. 

Notice, finally, that in this example the non-universality of the Higgs masses, given 
by (3.16), was chosen to be the maximal allowed by the modular weights {uh^ = — 1 
and hhj^ = —3). Also, the stau mass, for which we have nL,e = —1, cannot be further 
increased and therefore the behaviour under the UFB-3 constraint cannot be improved. 
Consequently, this optimised scenario represents a good estimate of how large the 
neutralino detection cross section can be in heterotic orbifolds with overall modulus, 
where soft masses are given by (3.11). We therefore conclude that in this class of 
models cr^o_p < 10~^ pb. The neutralino in these scenarios would escape detection in 
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Figure 10: The same as Fig. 6 but for the optimised example where = n^j^ = —1 with 
tan/? = 20, 35. 

all present experiments and only 1 tonne detectors would be able to explore some small 
areas of the allowed parameter space. 

For completeness, we have also analysed other three scenarios, described in [50], 
which also give rise to gauge coupling unification with an overall modulus. Their 
corresponding modular weights are summarised in Table 1. 

For example, in case C) unification is possible with ReT = 7, but extra massless 
chiral fields (one octet, one triplet, and two multiplcts transforming like right-handed 
electrons), with modular weight equal to —1, are needed. This scenario seems promis- 
ing, since the modular weights for sleptons are less negative {ni^eii ~ fact, 
although squarks become tachyonic at the GUT scale for cos^^ > 1/2, sleptons have 
a positive mass squared in the whole remaining area 6 G [vr/4, 3n/4], [57r/4, Tn/A]. As 
we have learned from the optimised example, this might be helpful in order to avoid 
the UFB constraints. The presence of extra matter alters the running of the gauge 
coupling constants, which are now dictated by the following beta functions, bi = —13, 
62 = —3, and 63 = 0. As a consequence, the running of the soft masses is also modified. 
In particular, all the gaugino masses become smaller at the EW scale, as compared to 
the usual running within the MSSM. Notice in particular that the gluino mass does 
not run (at tree level) from the GUT to the EW scales. 

The decrease in Mi and (especially) M2 affect the running of the scalar mass pa- 
rameters, rendering them smaller at the EW scale. This is enough to offset the increase 
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Figure 11: The same as Fig. 2 but for case C) in Table 1. 
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Figure 12: The same as Fig. 4 but for case C) in Table 1 and for tan/3 = 10 and 35. 
Notice that for tan/5 = 10 the lightest stop is the LSP on the dark grey area whereas for 
tan/3 = 35 it is the lightest stau. 

in mj^ ^ due to the smaller modular weights. Similarly, the important decrease in the 
gluino mass implies a very light squark sector. This leads to a qualitatively different 
structure of the SUSY spectrum in which squarks and sleptons have a similar mass. 
This is illustrated in Fig. 11, where we have represented the resulting spectrum for 
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Figure 13: The same as Fig. 5 but for case C) in Table 1 and with tan P — 10 , 35. 

= 200 GeV and tan/3 = 10 and 35. Unlike the previous cases, this example 
displays very light gluinos and squarks. There are even regions where the stop is the 
LSP (especially for low values of tan /3 for which the top Yukawa is larger) . 

The regions allowed by experimental constraints become larger in this example, as 
we can see in Fig. 12, where the {7113/2, 0) plane is depicted for tan/3 = 10 and 35. It 
is important to mention that, due to the resulting light squarks, the supersymmetric 
contribution to B{b — > 57) becomes sizable. Unlike in the previous examples, the 
experimental bound on this observable becomes the most stringent constraint, even for 
small values of tan /3. There are also areas which reproduce the correct dark matter 
relic density through coannihilation effects with the stop (for small values of tan/3) 
and the stau (for tan /3 > 20). Noticeably, in spite of the less negative modular weights 
for sleptons, the modifications in the RGEs (especially the decrease in Mi^2) imply 
smaller values for the Higgs mass parameters. In particular, mjj^ is more negative, 
making it more difficult to avoid the UFB-3 constraint. Only for heavy gravitinos do 
allowed regions occur (777,3/2 > 650 GeV is necessary for tan/3 = 10, whether for larger 
tan /3 gravitinos heavier than 1 TeV are needed) . As already observed in case A) , these 
regions never correspond to those with the correct neutralino relic density. 

Finally, it is worth mentioning that, since the Higgs mass parameters have a smaller 
departure from universality, we do not expect large neutralino detection cross sections 
in this example. The results for cr^o_p are represented in Fig. 13 for tan/3 = 10 and 35 
and, clearly lie beyond the reach of current and projected direct dark matter searches. 
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Figure 14: The same as Fig. 2 but for case D) in Table 1. 



Another potentially interesting scenario is case D), once more due to the reduced 
modular weights for sleptons. As in the previous example, the region allowed at the 
GUT scale is ^ e [n/i, Sn/i], [57r/4, 7n/4], where mQ^ j^^ > 0. In this scenario one 
needs four extra multiplets, transforming like {Q, Q, D, D), and Re T = 9 for unification 
to take place, thus implying bi = —12, 62 — —4, and 63 = 0. The absence of running 
for as and therefore for the gluino mass parameters has the same consequences as 
in case C), leading to a light squark spectrum. This is illustrated in Fig. 14, where 
the sparticle masses arc plotted as a function of the Goldstino angle for 7723/2 = 200 
GeV and tan/? = 10 and 35. Interestingly, the running of the wino mass parameter 
is slightly enhanced and Mi ~ M2 is found at the elcctrowcak scale. This makes the 
lightest ncutralino a mixed bino-wino state and almost degenerate in mass with the 
lightest chargino. Although this can lead to an increase of the resulting ncutralino 
direct detection cross section, it also implies a more efficient ncutralino annihilation 
and, consequently, a relic density which is too small to account for the dark matter 
of the Universe. In particular, one obtains VL^oH^ < 0.01 for the whole region with 
gravitinos fighter than 1 TeV, independently of the value of tan p. Thus, although 
the area allowed by experimental constraint, represented in Fig. 15 for tan (3 — 10 and 
35 is sizable, the astrophysical constraint on the relic density is never fulfilled. As 
in the previous examples, the presence of light squarks induce larger contributions to 
B{b — > 57) and the experimental constraint on it excludes extensive regions of the 
parameter space, even at low tan /3. Furthermore, the steeper running of M2 renders 
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Figure 15: The same as Fig. 3 but for case D) in Table 1. Notice that for tan/3 = 10 the 
lightest stop is the LSP on the dark grey area whereas for tan /3 = 35 it is the lightest stau. 

nijf^ more negative, which makes the UFB constraints even more restrictive. All the 
area represented in Fig. 15 becomes excluded for this reason. 

Finally, case E) in Table 1 corresponds to a Z'^ orbifold with a universal modulus, 
in which case ReT ~ 24 is needed. Due to the small modular weights for sleptons this 
example yields similar results regarding the UFB constraints as case A) , with allowed 
regions appearing only for very massive gravitinos and incompatible with the astro- 
physical bound on the dark matter relic density. Moreover, since the non-universality 
in the Higgs mass parameters {rrij^^ < m^J is not the optimal to increase the neu- 
tralino detection cross section, the theoretical predictions for cr^o_p are even smaller 
than those represented in Fig. 5. 

The presence of different moduli can provide some extra freedom to the non- 
universalities of the soft scalar masses (3.9), which are then parametrized by new 
Goldstino angles, ©j. We have analysed two scenarios of this kind, whose modular 
weights are described in Table 2, and which have also been shown to reproduce gauge 
coupling unification [50] . Scenario F) corresponds to a Zq orbifold with a rotated plane 
for which unification is achieved with ReTi ^ 10 ReT2. In scenario G), a, Z2 x Z2 
orbifold was taken, again non- isotropic, with Ti ^ T2,3. 

In both examples, some combinations of the Goldstino angles lead to potentially 
interesting non-universalities in the soft scalar masses. In particular, it is always pos- 
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Table 2: Modular weights for the scalar fields of two multimoduli heterotic orbifold scenarios 
that can reproduce gauge unification [50]. 

sible to enlarge the stau mass, thus avoiding UFB constraints. For instance, in case 
F) this can be done by choosing Q\ = 1 (hence Q\ = Q\ = 0) obtaining the following 
expressions for the soft terms. 
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(3.21) 


G) one can take 


= (and therefore Oi + e| = 


1) and obtain 




= ^3/2 (1 - ^ ®2 cos^ ^) , 














(3.22) 



Notice, however, that the resulting Higgs mass parameters are not adequate to obtain 
large neutralino detection cross sections. On the one hand, in case F) these are universal 
by construction since they have the same modular weights. On the other hand, in case 
G),the Higgs soft masses are xelaXcdhy m^^/m^j^ = (1 — 362 cos^ 6^)7(1 — 3 6^ cos^6'). 
If heavy staus are required, as in (3.22), m|^^ > Tnl^ is obtained, which implies small 
(TjjO-p- In order to reproduce the optimal Higgs non-universality, m|^^ > Tn\^, one 
needs to take 0^ > 0^ but then light staus appear (since m^^ = m|,^ = "^l^d) which 
therefore lead to strong UFB constraints. 

As a consequence, none of these two examples can provide larger theoretical pre- 
dictions for than those obtained in the optimised case of Fig. 9. 
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3.1 Non-thermal production of neutralinos from late gravitino 
decays 



Due to their extremely weak (gravitational) interactions, gravitinos can have a very 
long lifetime (typically longer than 10^ sec for 777.3/2 < ITeV). Consequently, gravitinos 
which are produced thermally in the reheating phase after inflation and which decouple 
with a given relic density 1^3/2 /^^ will decay at late times into the LSP (or cascading 
down to the LSP). This constitutes a source of non-thermal production of neutralinos, 
Q^o^^, which contributes to the total relic density of dark matter, and is related to the 
gravitino relic density as 

n^.rh' = ^ . (3.23) 

^1 7773/2 



It is easy, however to argue that Q^o is typically very small in the orbifold models 
we have just presented. It can be shown that in the case of gravitinos that are much 
lighter than gluinos (7773/2 <^ 777^), which generally holds in our case, the gravitino 
abundance from thermal production can be estimated as [67] 

^ ,2 f Tr \ /100GeV\ / 7715 \2 ^ 

The reheating temperature, Tr, is here a free parameter. The larger the reheating 
temperature, the more sizable $13/2/7^ becomes. However, in order not to reintroduce 
the infamous gravitino problem and spoil the results from Big Bang nucleosynthesis, 
the reheating temperature is constrained to be [68] < 10® GeV for a gravitino with 
a mass of order 7773/2 ~ 100 GeV. This implies an upper limit on the gravitino relic 
density, and thus on the non-thermal production of neutralinos. From (3.23) and (3.24), 
and using typical values for the gluino and neutralino mass, we obtain ^^(f^ ^ 2 x 10~^. 
Thus, even if the reheating temperature was increased to Tr = 10^ GeV (for heavier 
gravitinos, and provided their hadronic branching ratio is small), the contribution to 
the neutralino relic density is negligible. 
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4 Modifications due to the presence of an anoma- 
lous U{1) 



As discussed in the Introduction, an anomalous U{1) [28] is usually present after com- 
pactification^. For example, it was found in [71] that only 192 different three-generation 
models containing the SU{3) x SU{2) x U{1)"' gauge group can be constructed within 
the orbifold with two Wilson lines. The matter content of 175 of them was analysed 
in detail and only 7 of them turn out not to have an anomalous U{1) associated. Let 
us recall that the presence of the anomalous U{1) is crucial for model building [29]. It 
generates a FI contribution to the D-term [28], breaking extra U{1) symmetries, and 
allowing the construction of realistic standard-like models [30-32,35-38]. Let us now 
discuss the possible contributions to soft scalar masses generated because of the FI 
breaking, and their effect on the neutralino direct detection cross section. 

The FI breaking induces additional terms to soft SUSY-breaking scalar masses^ 
due to F-terms, namely, the so-called D-term contribution [39-44]. In particular, the 
presence of an anomalous U{1) after compactification generates the dilaton-dependent 
FI term. That is, the D-term of the anomalous U{1) is given by 

= ^fg:; + + TT'-qM , (4.25) 



where the first term corresponds to the dilaton-dcpcndcnt Fayet-Iliopoulos term with 
'gs 



the GS coefficient 6qq proportional to the value of the anomaly, and the second one is 



the usual D-term with the U{1) charges of the fields Then, some of these fields 
(with vanishing hypcrcharges), let us call them Cfs, develop large vacuum expectation 
values (VEVs) along the D-fiat direction in order to cancel the FI term, inducing the 
D-term contribution to the soft scalar masses of the observable fields. The resulting 
scalar mass squared is given by [41]^ 

-a=-3/2|l + -aCOS ^ + Z,{T + T*M<lirW J' ^ ^ ^ 



^In [70], some conditions for the absence of the anomalous U{1) are discussed, and classifications 
of models with anomalous U{1) are also attempted. 

^Let us remark that there are no additional contributions to gaugino masses and A-terms when 
Higgs fields relevant to such symmetry breaking have less F-term than those of dilaton and moduli 
fields. 

^Assuming for simplicity that the fields C/3 have no other U{1) charges. 
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A-I) 


-2 


A-II) 


-2 -2 


A-III) 


-2 1 -2 



Table 3: Variations of scenario A) in Table 1. The fields have the same modular weights, 
but we have assigned the above U{1) charges to the slepton and Higgs fields. 

where the first two terms are the usual contributions (see (3.11)), and the third term 
is the D-term contribution. Obviously, if the observable fields have vanishing U{1) 
charges, q^, this contribution is also vanishing, and we recover the situation of the 
previous Section. However, the observable fields have usually non- vanishing charges in 
explicit models [30-38] , and the effect of this contribution must be taken into account 
in the analysis. 

As wc can see in the above formula, the D-term contribution generates in general 
an additional non- universality among soft scalar masses which depends on q^. Let us 
simphfy the analysis considering the case where only a single field C develops a VEV 
in order to cancel the FI term. Thus the above result reduces to the following form: 

= ^3/2 1 1 + cos' ^ + ^ [(6 - nc) cos' ^ - 5] I , (4.27) 

where q^ and ric are the U{1) charge and modular weight of the field C, respectively. 
It is worth emphasizing here that even in the dilaton-dominated case (cos^ = 0) the 
soft scalar masses are non-universal, 

^'=^3/2(1-5^) • (4.28) 

This result should be compared with the one in (3.11). The contributions —^q^/q^ 
correspond to the 5's in (2.1) and (2.2). It is noteworthy that, contrary to the cases 
without an anomalous U(l) analysed in the previous Section, positive values for the 
non-universalities are now possible, by choosing qa/lc < 0- This is welcome, as we 
will see, in order to enhance the stau masses and thus avoid the UFB constraints, or 
increase the value of ra\^. Moreover, we may expect that, for appropriate values of the 
U{1) charges, the above additional terms lead to interesting values for the neutralino 
detection cross section. 

Let us consider again scenario A) of the previous section, defined by expressions 
(3.15) and (3.16), where the sleptons masses were too small and therefore problematic. 
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Figure 16: The same as Fig. 2 but for example A-l) of Table. 3 with tan/3 = 35 and 45. 
Soft masses satisfy (3.15) and (3.16) with slepton masses modified according to Eq.(4.30). 



Using (4.27), these are now given by 



ml„e, = ml/, 1 1 - 3 cos^ 9 + ^ [(6 - nc) cos' ^ - 5] | . (4.29) 

For instance, assuming that C is a twisted field with modular weight nc — —2, and 
the ratios^° Q.Li^,eJ^c — ~2, one has 

^4,e« = ^3/2(ll-19cos^^) . (4.30) 



This corresponds to case A-I) in Table 3. The degree of non-universality, as defined in 
(2.2), is therefore (^L^.e^ = 10 — 19cos^ 9. Unlike the examples in the previous section, 
where all the 5's had to be negative, positive values can now be obtained. 

Condition cos^ ^ < 1/3 is no longer necessary in order to avoid tachyonic sleptons, 
since cos^ ^ < 11/19 is sufficient, thus implying broader allowed regions. Noticeably, 
unlike the examples in the previous Section, gaugino masses (3.17) can be of the same 
order, and even smaller than slepton masses at the GUT scale. For example, with 
cos^^ = 1/2 one obtains m\^,,^ — with gaugino masses of the same order 



^°These values for modular weights, and ratios between ?7(1) charges of observable and FI fields, 
are typically obtained in explicit models [30-32] . 
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~ mj ^ . And in the dilaton dominated limit one obtains m? „ =11 7n^,r,, much 
larger than gaugino masses, for which ~ 37713^2- 

Very heavy sleptons are therefore possible and, although their masses still decrease 
away from the dilaton limit, in most of the parameter space their mass is similar to 
that of squarks. Such a heavy scalar sector is clearly shown in Fig. 16 for 7/13/2 = 200 
GeV and tan/? = 35 and 45. As we can sec, the lightest neutralino and chargino 
are much lighter and the lightest stau rarely becomes the LSP. Only along extremely 
narrow areas, which appear for tan /5 > 25 at the edge of the allowed zones, is the 
stau lighter than the neutralino. Consequently, the allowed parameter space where the 
lightest neutralino is the LSP is more extensive and larger values of tan j3 are allowed. 
This leads to a more efficient decrease of both the pseudoscalar and heavy scalar Higgs 
masses. 

The resulting (7713/2, 0) parameter space is represented in Fig. 17 for tan = 35 , 45. 
The most salient consequence of the increase in the slepton masses is the reduction of 
the regions excluded by the UFB-3 constraint (this is particularly true in the dilaton 
limit, as pointed out by [72]). Contrary to what happens in case A), for which the whole 
parameter space was excluded (see Fig. 4), now only some points with 7/23/2 < 100 GeV 
are disfavoured. On the other hand, as a consequence of the increase in tan/5 and 
the associated decrease of the pseudoscalar Higgs mass, the experimental bounds on 
H'^ and b ^ s'y become much more constraining. These imply a stronger 
lower limit on the gravitino mass as we can see in the plot. 

Regarding the neutralino relic density, for tan (3 < 25, the stau coannihilation re- 
gions are not present and ^j^oh? > 0.3 is obtained in the whole (7/23/2, 0) plane. The 
coannihilation strip only appears for tan /? > 25, where the correct relic density can 
be reproduced in extremely narrow areas. However, the most interesting case is 
tan f3 > 30, where thanks to the decrease in the pseudoscalar Higgs mass, we can find 
regions in the parameter space where ^m^o ^ 171^, in which case the resonant anni- 
hilation of neutralinos helps reproducing the WMAP results. Due to this effect, the 
regions of the parameter space satisfying experimental and astrophysical constraints 
are now significantly larger. For instance, with tan (3 — 35 we find allowed points 
for the whole range of gravitino masses between 350 GeV to 1000 GeV. Notice also 
that in the dilaton hmit, 9 — 7r/2, the region reproducing the WMAP rehc density is 
experimentally excluded for tan P < 47. 

The decrease on the Higgs masses for tan (3 >35 leads to an enhancement of the 
neutralino-nucleon cross section. However, for such values of tan/3 large branching 
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Figure 17: The same as Fig. 3 but for example A-l) of Table. 3 with tan /3 = 35 , 45. Soft 
masses satisfy (3.15) and (3.16) with slepton masses modified according to Eq.(4.30). 
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Figure 18: The same as Fig. 5 but for example A-l) of Table. 3 with tan /3 = 35 , 45. Soft 
masses satisfy (3.15) and (3.16) with slepton masses modified according to Eq.(4.30). 



ratios for the rare processes 6 — > 57 and — > /x"*"//" are found, which therefore put 
stringent upper bounds on cr^o_p. This can be seen in Fig. 18, where the theoreti- 
cal predictions for the neutrahno detection cross section are represented as a func- 
tion of the neutrahno mass for tan/3 = 35, 45. The cross section is bounded to be 
(7^0 _p < 2 X 10~^ pb by the aforementioned constraints. The corresponding values for 
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Figure 19: The same as in Fig. 6 but for example A-l) of Table. 3 with tan = 35 , 45. Soft 
masses satisfy (3.15) and (3.16) with slepton masses modified according to Eq.(4.30). 

H and pseudoscalar mass are shown in Fig. 19. Both are rather large, comparable to 
what we obtained in Fig. 10 for the optimised case in the previous Section. For such 
large tan /3 the small CP-odd Higgs mass often leads to an excessive contribution to 
Bg — > n'^n", incompatible with the present experimental bound. 

In order to further increase the neutralino-nucleon cross section while avoiding the 
experimental constraints, we need to find a way to decrease the Higgs masses even for 
small values of tan /3. This can be done by further exploiting the D-term contribution 
if we assume that also the Higgses have non- vanishing anomalous U{1) charge. For 
example, setting q^^/qQ = —2, using result (4.27) equation (3.16) is modified as 

= "^3/2 (11 - 17cos2 e) , 

= (1-3 cos^ 0) . (4.31) 

We have labelled this as case A-II) in Table 3. The degree of non-universality, using 
notation (2.1), is given by Sh^ = 10 — 17cos^^^ and = — 3cos^^. In the dilaton- 
dominated case this turns out to be very large, Sh^, = 10. Notice once more that, due 
to the presence of the anomalous U(l), positive values for 5h„ can be obtained. 

Due to the large non-universalities, the fi parameter and heavy Higgs masses can be 
significantly reduced, even for moderate values of tan/?, and especially in the dilaton 
limit. This is evidenced in Fig. 20 where the SUSY spectrum is plotted for 1713/2 = 
200 GeV and tan/3 = 10 and 25. The lightest neutralino has an increased Higgsino 
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Figure 20: The same as Fig. 2 but for example A-ll) of Table. 3 with tan/3 = 10 and 25. 
Soft masses satisfy (3.15) with slepton masses modified according to Eq.(4.30) and Higgs 
mass parameters given by (4.31). 

composition in the area close to the dilaton-dominated hmit and its mass is similar to 
that of the lightest chargino. This, together with the decrease in m^o, leads to a more 
effective ncutralino annihilation in the early Universe and hence a reduced relic density 
towards — 7r/2. 

The associated 0) plane is represented in Fig. 21 for tan/3 = 10 and 25. 

Extensive areas of the parameter space become excluded since m\ becomes negative. 
These regions correspond to the ruled areas above the experimentally allowed ones. 
They are specially constraining in the dilaton limit, for which m3/2 ^ 300 GeV is 
needed. As tan j3 increases the dilaton limit eventually becomes excluded and the 
allowed regions shrink towards 6 = 7r/4 and Sir/ 4. As mentioned above, the resulting 
neutralino relic density is too small in the dilaton limit, but away from it, regions 
appear where the WMAP result is reproduced. Once more, due to the increase of the 
stau mass and rrijj^, the UFB constraints pose no problem in this scenario, and all the 
points depicted in Fig. 21 satisfy them. 

The decrease of the parameter and the resulting larger Higgsino component for the 
lightest neutralino, A^^3+A^^4 ^ 0.4, contribute to the increase of o'^o_p. The neutralino- 
nucleon cross section is depicted on the left-hand side of Fig. 22 for tan f3 = 25, display- 
ing neutralinos with (7^o_„ > 10~^ pb and a mass in the range m^o 200 — 500 GeV 
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Figure 21: The same as Fig. 3 but for example A-ll) of Table. 3 with tan/3 = 10 and 25. 
Soft masses satisfy (3.15) with slepton masses modified according to Eq.(4.30) and Higgs 
mass parameters given by (4.31). 
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Figure 22: The same as Fig. 5 but for example A-ll) of Table. 3 with tan/3 = 25. Soft 
masses satisfy (3.15) with slepton masses modified according to Eq.(4.30) and Higgs mass 
parameters given by (4.31). 



that would be within the reach of the CDMS Soudan or XENONIO experiments. On 
the right-hand side of Fig. 22 the corresponding value of the // parameter is represented 
as a function of the pseudoscalar Higgs mass for the same example, evidencing the de- 
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crease in both quantities. The points with a larger detection cross section correspond 
to those with ttia ^ 220 - 240 GeV and f« 350 - 500 GeV. 

Notice that, unlike the previous case, a significant increase of cr^o_p is obtained 
while keeping the observables B{b — 57) and B(i?^ /^-^/^ ) under control. As 
already mentioned, this owes to the fact that large values of tan /3 are no longer needed 
for obtaining sizable detection cross sections. 

The Higgs non-universality can be further increased if we also consider to be 
charged under the anomalous U{1). For example, taking qjj^/qc = 1/2, which corre- 
sponds to case A-III) in Table 3, the soft masses for the Higgses become 

= m^/2(ll-17cos2e) , 

ml, = ml/2 (-^ + cos^9^ . (4.32) 

In this case the non-universalities are given by 5^^^ = 10 — 17cos^^, and Sh, — — § + 
cos^ 9, which implies that in the dilaton limit they become 5h^ — 10 and 5h, — ~|- 

The associated parameter space is represented in Fig. 23 for tan/3 = 20. As in 
the previous example, extensive areas of the (777.3/2, 0) plane are ruled out because of 
777^ becoming negative. Now, the region around the dilaton-dominated case, where 
the non-universality is maximal, is completely excluded for this reason. Only narrow 
allowed areas of the parameter space for specific values of the Goldstino angle survive. 
These allowed areas become smaller as tan/3 increases, since the pseudoscalar mass 
becomes tachyonic more easily, and eventually disappear for tan /3 > 30. 

As a consequence of the further decrease in 777^^^, even smaller values of tan/3 are 
enough in order to obtain light pseudoscalar and scalar Higgses. Consequently, values 
of cr^o_p within the reach of present experiments can now be found for even smaller 
values of tan/3. These are shown as a function of the neutralino mass on the left- 
hand side of Fig. 24 for tan /3 = 20, which is the optimal choice for this case. The 
pseudoscalar mass is very small in these points, as evidenced on the right-hand side of 
Fig. 24, where the /i parameter is plotted versus 777^- The points with (T^o_p > 10~^ pb 
correspond to those with 777^ < 200 GeV, and have a moderate Higgsino composition. 

Larger values of qH^/lc would further reduce 777^^^ and lead to more sizable non- 
universalities. However, this makes it more difficult to obtain positive fi^ or 777^, and 
the allowed areas in the parameter space become more fine-tuned. For instance, with 
lHa/lc ~ 1' which leads to 

"^k = "^3/2 (11 - 17cos^ e) , 
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Figure 23: The same as Fig. 3 but for example A-lll) of Table. 3 with tan/3 = 10 and 20. 
Soft masses satisfy (3.15) with slepton masses modified according to Eq.(4.30) and Higgs 
mass parameters given by (4.32). 
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Figure 24: The same as Fig. 5 but for example A-lll) of Table. 3 with tan/3 = 20. Soft 
masses satisfy (3.15) with slepton masses modified according to Eq.(4.30) and Higgs mass 
parameters given by (4.32). 



rriu^ = mli2 (-4 + 5 cos d) , 



(4.33) 



there are only two narrow regions left with 200 GeV < 171^/2 < 400 GeV around 6 ^ 0.8 
and 9 pa 2.3, for which ^L^.e^ ~ 0.8, 6h^ 1.7, 6hj^ ~ -2.6, and ^L^.e^ ~ 1-6, Sh^, ~ 2.5, 
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Figure 25: The same as Fig. 5 but for tan/? = 10, 15, 20 and 25. Soft masses satisfy (3.15) 
and the effect of D-terms on the slepton and Higgs masses has been included according to 
Eqs.(4.3G) and (4.33), respectively. 



5hj^ ~ —2.8, respectively. The results for cr^o_p are qualitatively similar to those of 
Fig. 24, with points within the reach of CDMS Soudan and XENONIO for tan /3 = 10 
to 25. The complete (m3/2, 9) plane is ruled out for tan/3 > 25 due to the smallness of 
m\. 

In order to illustrate this possibility, the theoretical predictions for o-^o-p arc repre- 
sented in Fig. 25 for a scan of the cases with tan j3 = 10, 15, 20 and 25. The neutralino 
detection cross section reaches the sensitivity of the CDMS Soudan and XENONIO 
experiments for m^o fa 150 — 200 GeV. On the right-hand side of Fig. 25 the corre- 
sponding values of are plotted as a function of the CP-odd Higgs mass, evidencing 
the important decrease in both as a consequence of Higgs non-universalities. 

If the ratio of f/(l) charges is further increased to Qh^/Qc = 2, the whole parameter 
space becomes excluded. 

Notice finally that, although we have exemplified the relevance of the contribution 
due to an anomalous U{1) with a variation of case A) of Table 1, the rest of the 
examples in that table can also improve their behaviour under UFB constraints and 
lead to detectable neutralino dark matter by using this procedure. This owes to the fact 
that the non-universality introduced by the D-term contribution typically dominates 
over the modular weights, especially in the dilaton limit, as we noted in Eq. (4.28). 



40 



5 Conclusions 



We have studied the phenomenology of supergravity theories which resuh from orbifold 
scenarios of the heterotic superstring and the resulting theoretical predictions for the 
direct detection of neutrahno dark matter. These scenarios are specially interesting, 
since the soft terms can be computed explicitly and have, in general, a non-universal 
structure. 

We have studied the parameter space of these constructions, computing the low- 
energy spectrum and taking into account the most recent experimental and astro- 
physical constraints. In addition, we have imposed the absence of dangerous charge 
and colour breaking minima. In the remaining allowed regions the spin-independent 
part of the neutralino-proton cross section has been calculated and compared with the 
sensitivity of current and projected dark matter detectors. 

In the absence of an anomalous U{1) the non-universality of the soft scalar mass 
parameters is always negative with respect to the gravitino mass. The smallness of 
the stau mass implies that the UFB-3 direction in the parameter space becomes very 
deep and the realistic minimum is no longer the global one. In most of the cases it 
is not possible to satisfy this constraint in the regions permitted by astrophysical and 
experimental bounds. Also, the non-universality on the Higgs masses is not sufficient to 
produce any significant increase of the neutralino-nucleon cross section. Consequently, 
the theoretical predictions for (J^o_p are typically beyond the sensitivity of dark matter 
experiments. 

The presence of an anomalous U{1) gives more flexibility to the non-universality 
in the scalar masses. These can even be heavier than the gravitino mass at the GUT 
scale. This allows to increase the slepton masses, thus avoiding the UFB constraints. 
Moreover, the lightest ncutralino becomes the LSP in most of the parameter space 
and regions fulfilling all the experimental and astrophysical constraints can be found. 
Furthermore, the non-universahty of the Higgs mass parameters can be tailored to 
favour the presence of fight Higgses and the increase in the Higgsino composition of 
the fightest neutrahno. This leads to a sizable increase in the theoretical predictions 
of cr^o_p and compatibility with present experiments can be achieved. 
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